Tunable Magnetic Relaxation In Magnetic Nanoparticles 
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We investigate the magnetization dynamics of a conducting magnetic nanoparticle weakly coupled 
to source and drain electrodes, under the assumption that all relaxation comes from exchange of 
electrons with the electrodes. The magnetization dynamics is characterized by a relaxation time 
ii, which strongly depends on temperature, bias voltage, and gate voltage. While a direct measure 
of a nanoparticle magnetization might be difficult, we find that t\ can be determined through a 
time resolved transport measurement. For a suitable choice of gate voltage and bias voltage, the 
magnetization performs a bias-driven Brownian motion regardless of the presence of anisotropy. 
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Electrical and magnetic dynamics of itinerant ferro- 
magnets are often described using a generalization of 
"Born-Oppenheimer" approximation, in which the elec- 
trons adjust instantaneously to any changes in the mag- 
netization. This approximation, which is justified by the 
wide disparity of electronic and magnetic time scales, is 
found to hold down to the smallest size scales attainable 
by present-day nanofabrication techniques. The list of 
physical phenomena it explains includes the Giant Mag- 
neto Resistance effect the exchange interaction be- 
tween different layers in ferromagnetic multilayers 0, 
and the nonequilibrium spin torque 0- The first-order 
effect of a changing magnetization on the electrons was 
considered only recently . 

Once electronic and magnetic degrees of freedom are 
separated, a realistic description of magnetization dy- 
namics requires the inclusion of a relaxation mechanism. 
A phenomcnological description of magnetic relaxation 
is provided by the Gilbert damping term in the Landau- 
Lifschitz-Gilbert equation. The Gilbert damping term 
represents the magnetic relaxation caused by spin-orbit 
scattering, phonons, etc; its microscopic origin is the sub- 
ject of ongoing research Additional magnetic relax- 
ation in thin ferromagnetic films follows from the emis- 
sion of spin currents into normal metals adjacent to the 
ferromagnet 0. 

It is the purpose of this letter to study time-dependent 
electric and magnetic properties of a magnetic nanoparti- 
cle, weakly coupled to source and drain reservoirs via tun- 
neling contacts. Such magnetic nanoparticles have been 
fabricated and studied recently by Ralph and cowork- 
ers [1, 0- In the absence of strong spin-orbit coupling 
in the nanoparticle or for sufficiently large conductances 
of the tunneling barriers, all electronic and magnetic re- 
laxation occurs via the exchange of electrons with the 
leads. This system is sufficiently simple that magnetic 
and charge degrees of freedom can be treated on equal 
footing and the separation of time scales for charge and 
magnetization dynamics can be derived from a micro- 
scopic model. Moreover, because magnetic relaxation 



takes place through the exchange of electrons with the 
leads, the magnetic dynamics crucially depends on the 
electric environment (bias voltage and the voltages on 
nearby metal gates). This dependence leads to a tunable 
magnetization relaxation rate. To our knowledge, this is 
the first system in which magnetic damping rate can be 
tuned by a simple gate voltage. Tunability is an impor- 
tant asset for potential applications of magnetic nanos- 
tructures, since optimal functioning of a nanomagnetic 
device requires that damping rate is matched to other 
relevant time scales of the system. Further, as we show 
below, for certain values of bias and gate voltages, the 
electric current creates a non-equilibrium "randomiza- 
tion" process that exceeds relaxation, causing a random 
motion of the magnetization vector, even in the presence 
of magnetic anisotropy. 

Model. Starting point of our discussion is a model 
Hamiltonian for a ferromagnetic nanoparticle, in which 
the ferromagnetism arises from the long-range exchange 
interaction 

H dot ^Y. e A^- J S-S-^S 2 z +E c {N-N g ) 2 . (1) 
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Here c^ a and cL, are annihilation and creation operators 
for an electron with spin a and energy e^, 
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is the total spin of the nanoparticle, a being the vector 
of Pauli matrices, J and K set the strength of the ex- 
change interaction and anisotropy, respectively, Eq is the 
charging energy, and N g is proportional to the gate volt- 
age. Without the anisotropy term, this model is known 
as the "universal Hamiltonian" , which has been shown 
to describe non-magnetic metal nanoparticles on energy 
scales below the Thouless energy (Note, however, 

that, unlike the "universal Hamiltonian", Eq. should 
be considered a model description only, since ferromag- 
netism implies a splitting between majority and minority 
electrons that exceeds the Thouless energy.) 
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FIG. 1: (a) Schematic of the system under consideration: 
a ferromagnetic nanoparticle (F) is connected via tunneling 
junctions to two electrodes, which can be normal (N) or fer- 
romagnetic (F) metals. A bias voltage Vbias is applied across 
the nanoparticle, while a gate voltage V g can be applied on 
a gate capacitively coupled to the particle, (b) Structure of 
the electronic ground state of the nanoparticle. The ground 
state has N a singly occupied levels and Fermi levels Em and 
Em for majority and minority electrons, respectively. 

Without exchange interaction and magnetic 
anisotropy, Eq. |QJ is the basis for the "constant 
interaction model" of the Coulomb blockade [10|, For 
sufficiently large exchange interaction, the ground state 
of the Hamiltonian (JIJ is ferromagnetic; it has N s = 2S 



singly occupied single-electron levels and different Fermi 
levels Em and E m for majority and minority electrons, 
respectively, see Fig. with Em — E m = 25 J. The 
excited many-electron states \a) are characterized 
by occupation numbers of the single-electron levels 
= 0,1,2, the total spin S, the z component of the 
total spin, 5 Z , and an additional degeneracy parameter 
if 5 is not maximal. (The z axis is assumed to be the 
easy axis of the magnet.) For low-lying excited states, 
the identity N s — 25 is preserved, since many-electron 
states with Na > 25 have an excitation energy of order 
Em — E m [ll|, |l3] • The typical energy scale for low- lying 
excited states is of the order of the exchange energy J. 
In the nanoparticles studied in Refs.0,0, the anisotropy 
K is of order J as well ^ij, so that the magnetic 
excitations S x — * S z ± 1 have energy comparable to that 
of electron-hole excitations. 

The time-independent low energy properties of this 
model have already been extensively analyzed in Refs. 
Iffl fl2l Here we study the time dependence of the total 
spin 5 for the case when the nanoparticle is connected 
to source and drain reservoirs, which can be normal met- 
als or ferromagnets, see Fig. QJi. If the leads are ferro- 
magnetic, we assume that the magnezation is collinear 
with the easy axis of the nanoparticle. For weak cou- 
pling (level broadening much smaller than temperature 
and level spacing), the electronic state of the nanoparti- 
cle can be characterized by the probability P a to find the 
nanoparticle in many-electron state \a). The probabili- 
ties are governed by a rate equation, 



dP a 
dt 



EE E E {r^\{a%>)\ 2 [(l-ME a -E«'))Pa-fi(E a -E a ,)P a ,] 
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Here, the sum over a' extends over all many-electron 
states \a'), whereas the sum over /i extends over all 
single-electron states [i. The distribution functions for 
the left and right reservoirs are denoted /;, I = L,R, 
whereas r ' p\ a set the tunneling rate through contact I 
for majority (|) or minority (!) electrons to the one- 
body state /1. Equation J3J generalizes the rate equation 
used in the theory of Coulomb blockade However, 
one should note that upon using an evolution equation 
for probabilities, information on off-diagonal terms of the 
density matrix is lost. This means that, while we keep 
track of the magnetization along the z-axis, the position 
of the magnetization in the xy plane is not monitored. 
This prevents the use of this approach to the case of fer- 
romagnetic leads with non-collinear magnetizations. 



Separation of time scales. The magnetic structure 
of the rate equation is encoded in the matrix elements 
(a'lc^la). Only matrix elements for which the particle 
number N a r — N a ± f and the spin Sq,> — S a ±1/2 are 
nonzero. The spin-dependence of the corresponding tran- 
sition rate (matrix element squared) is proportional to 
the square of a Clebsch-Gordon coefficient Close to 
equilibrium (S z ~ S), these two relevant Clebsch-Gordon 
coefficients are of order 1 and 1/V~S <^ 1, for up and down 
electrons respectively tunneling on a majority state, and 
vice versa for a minority state. As we shall see below, 
it is this hierarchy that leads to the separation of time 
scales for the charge and magnetic degrees of freedom. 

We now consider Eq. (0 in the limit where both the 
temperature and bias voltage are much smaller than J, 
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so that only one single-electron state /j, is involved in 
transport, and for the case that the system is close to 
equilibrium, n = S—S z <C 5. For defmiteness, we assume 
that fi is a minority state. Denoting the many-electron 
state with N electrons, spin 5 and S z = S—n by "(n, — )" 
and the many-electron state with N + 1 electrons, spin 
5 - 1/2 and S z = S - 1/2 - n by "(n, +)" , Eq. © reads 



d,P+ = 



E 

l=L,R 



P H 



dtP- = 



l=L,R 



d: 



(4) 
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and we defined //„x = - E n ), and //„ T = //(£"+ - 

E~,i). The structure of the rate equation is similar if 
transport is facilitated by a majority state, or if several 
one-electron states contribute to transport. In all cases, 
one retains the special hierarchy of transition rates that 
is present in Eqs. (JSJ) and (JSJ: n = S — S z stays constant 
in most tunneling events (transition rate of order r„;jj, 
whereas only a small fraction of tunneling events allows n 
to change by unity (transition rate of order T^n/(2S + 
!))• 

To leading order in n/S, only the second lines in Eqs. 
(JSJ and ©, corresponding to a process where a down 
spin enters or exits the quantum dot, contributes. Such 
a process facilitates charge relaxation, the corresponding 
relaxation time scale being 



U = (r M ix + TpRi) 



(7) 



This is the same time scale as for charge relaxation in a 
non-magnetic Coulomb blockaded particle . 

The magnetization dynamics appear when we consider 
terms of order n/S in Eqs. JSJ and ©. We consider the 
n-dependence of the Clebsch-Gordon coefficients in Eqs. 
(J5J and |JBJ, but neglect the (very weak) n-dependence 
of the Fermi functions f, and f,„ . Assuming that the 
charge degree of freedom is in equilibrium, the probability 
for the nanoparticle to have spin S z = S — n is P n = 
P r f + P,7 = (1 + c)P~, where c is a weighted average of 
the distribution functions in the two leads, 



(8) 



(The ratio c = / P n is taken n- independent in view of 
the comments made above.) The evolution equation for 



P n then reads 

d t P n = 7-( n + + 7+nPn-i 

- [ 7 -n + 7+ (n+l)]P n , 

where we abbreviated 



(9) 



7+ 



7- = 



25+1 

(1 - c)T^fl + (1 - c)r MflT /T 
25 + 1 



(10) 



In equilibrium one has 7_ > 7+, corresponding to the 
stable stationary solution 



Pn 



( 7+ /7_)»(l-7 + / 7 _). 



(11) 



Equation J^J has a simple interpretation in terms of a 
random walker on a semi-infinite chain that has a prob- 
ability 7_rt to step towards the origin and a probabil- 
ity 7-|_ (n + 1) to step away from the origin. The dif- 
ference 7_ — 7+ can be identified as the magnetic re- 
laxation rate l/£i; the sum 7+ + 7_ corresponds to a 
"randomization rate" for the magnetization. The iden- 
tification 7- — 7+ = 1/ii becomes manifest when the 
time-dependent problem is cast in terms of an evolution 
equation for moments of P n (corresponding to moments 
of S z ). For the first moment A4(t) — X)n n -^"W one 
finds 



B t M(t)='Y + -M(Jt)/t 1 , t- 1 =7_- 7+ . 



(12) 



which is the equivalent of the Landau-Lifschitz-Gilbert 
equation for our system. Note that, while the charge 
relaxation rate 1 /t c depends on the tunneling rates only, 
both 7_ and 7+ depend strongly on temperature, bias 
voltage, and gate voltage via the Fermi functions /£ and 
fn- 

If a bias voltage is applied on the system, the effect of 
the relaxation of M. (t) toward its equilibrium value can 
be inferred from a small change of the measured electric 
current, thus allowing for an electric measurement of the 
magnetic relaxation time t\ , 



I = 



r 



eM(t) 

2S +1 { Tf+L i + T uR i 



t^Rl M-^T (j! 



Mil 



■(/i-c)- W/i-c) 



(13) 



[In Eq. (|13fl constant terms of order 1/(25+1) have been 
dropped.] Whereas additional intrinsic relaxation mecha- 
nisms (phonon-magnon coupling or spin-orbit coupling), 
which have not been included in the above analysis, may 
alter t\, Eq. I|13|l remains valid in all cases. 

Non equilibrium induced Browian motion of the mag- 
netization. An interesting limit occurs when when the 
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FIG. 2: Current / (upper) and reduced magnetization {S z )/S 
(lower) as a function of the bias voltage Vbiaa for a system 
with J = 1, K = 0.1 and k B T = 0.01. The down (up) 
arrows correspond to minority (majority) states that become 
involved in the tunneling processes. 

bias voltage is bigger than Eq — Eq , whereas the tem- 
perature is much smaller than Eq — Eq . In that case, 
one may set ff n = ff n = 1 while ff n = ff n = 0. For 
normal metal leads and a symmetric particle-lead cou- 
pling {T La — Tflcr), one then obtains that the station- 
ary solution is given by the detailed balance solution, 
= P~ = constant, independent of the magnetization 
n. [This result can be verified from direct solution of 
the detailed balance equation associated with Eq. ©]. 
This total randomization of the magnetization occurs ir- 
respectively of the presence of the anisotropy: although 
the (anisotropy) energy of a state increases with n, its 
probability P n does not decrease. This remarkable result 
is closely linked to the magnetic structure of the mas- 
ter equation: The bias voltage serves as a driving force 
for magnetization randomization, whereas the magneti- 
zation relaxation is suppressed by Coulomb blockade. We 
note that through this mechanism, the particle's mag- 
netic energy can increases to a value much above the 
ground state energy. In a real sample, this effect is lim- 
ited by extra sources of relaxation that we have not in- 
cluded (e.g., spin-orbit scattering). 

In Fig. |21 an example of the bias-voltage dependence 
of current and average magnetization is shown, obtained 
from numerical solution of the rate equation @ . Except 
small regions of width AVbi as = K/e near current steps, 
which appear when a transition n + 1 — > n is allowed 
while the corresponding n — > n transition is still blocked, 
one observes that the magnetization is fully randomized 
as soon as the bias voltage exceeds the threshold for cur- 
rent flow. The randomization effect is suppressed when 
K > J. At this point, particle- hole excitations affect the 
hierarchy in the structure of the rate equation @ , allow- 
ing for additional relaxation mechanisms and the break- 



down of detailed balance. (For example, an excited state 
with N + 1 electrons and n = 0, which can be accesses 
from the iV-electron ground state with n = 1, can relax 
to the iV-electron ground state with n = 0. However, the 
opposite process, a direct transition from the iV-electron 
ground state with n = to an excited N + 1-elcctron 
state with n = 1, is not allowed.) 

Conclusion. We have studied a simple but realistic 
model of a magnetic dynamics in a magnetic nanoparti- 
cle, for which relaxation occurs via the exchange of elec- 
trons with source and drain reservoirs. For this system, 
electric and magnetic properties can be studied on the 
same footing, and the separation of electric and magnetic 
time scales can be established explicitly. The magnetic 
dynamics is characterized by relaxation and randomiza- 
tion rates, which can both be tuned by bias voltage and 
gate voltage. A damping rate well matched to other time 
scales is crucial for magnetic dynamics that is both fast 
and robust. In this sense, the tunability of the damping 
rate discovered here may be of use in future studies of 
magnetic dynamics on the nanoscale. In addition, the 
possibility of measuring t\ through a transport measure- 
ment removes the need of a direct magnetization mea- 
surement, which is difficult for nanoscale magnets. 
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